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Abstract 

Travel-times of acoustic signals were measured between a hot tom- mounted source near 
Oahu and four bottom-mounted receivers located near Washington, Oregon, and Cali- 
fornia in 1988 and 1989. This paper discusses the observed tidal signals. At three out 
of four receivers, observed travel times at M2 and S2 periods agree with predictions 
from barotropic tide models to within ±30° in phase and a factor of 1.6 in amplitude. 
The discrepancy at the fourth receiver can be removed by including predicted effects of 
phase- locked baroclinic tides generated by seamounts. 

Our estimates of barotropic M2 tidal dissipation by seamounts vary between 2 x 10 16 
and 1 X 10 18 erg-s" 1 . The variation by two orders of magnitude is due to uncertainties in 
the numbers and sizes of seamounts. The larger dissipation (1 X 10 18 erg-s" 1 ) is the same 
order as previous estimates and amounts to 4 % of the total dissipation at M2. 

Thesis Supervisor: Dr. John L. Spiesberger 
Woods Hole Oceanographic Institution 
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Chapter 1 



Introduction 



We measure travel-times of acoustic signals between a bottom-mounted source near Oahu 
and four bottom-mounted receivers located near Washington, Oregon, and California 
in 1988 and 1989 (Figure 1-1). The power-spectra of basin-scale travel-times exhibit 
prominent variability at tidal periods [Spiesberger et al., 1989b]. There are two important 
reasons for further study of the tidal signals. First, numerical models for barotropic 
tides can be validated by using these new observations. Second, it is important to detide 
tomography data to better study other ocean fluctuations at similar periods [Spiesberger 
et al., 1989b]. 

The tides can be separated into barotropic and baroclinic modes, and their associated 
influences on travel-time can be separated into current and displacement effects. The tidal 
currents affect travel-time by alternating between alliance and opposition to the direction 
of sound propagation. The rise and fall of the water column affects the travel-time by 
causing periodic vertical iso-therm displacements and pressure fluctuations. Worcester 
et al. [1990] has shown that barotropic current effects can be quantified by using recip- 
rocal transmissions; one-way transmissions are used in this experiment, so we quantify 
barotropic effects by comparing measured data to barotropic tide model predictions. 

We find that acoustic travel-times computed from barotropic tidal models are similar 
to observations at three out of four receivers. The discrepancy at the fourth receiver 
can be understood by including the effects of phase-locked baroclinic tides generated by 
seamounts near the transmission path. Generation mechanisms are important both for 
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Figure 1-1: The acoustic tomography experiments of 1988 and 1989 consist of a source 
located near Kaneohe Bay, Oahu, and four bottom receivers (III to R4) whose approximate 
positions are shown. The distances between source and receivers are all about 4000 km. 
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understanding the distribution of internal wave energy and the energy budget of tidal 
dissipation [Munk, 1966; IVunsch, 1975; Bell, 1975; Lambeck, 1980; Hendershott, 1981]. 

This report closely follows Headrick et al. [1990] and is organized as follows. Chapter 
2 discusses the experiment and the estimation of acoustic travel-time. Comparisons be- 
tween observations and predictions from barotropic and baroclinic tide models are made 
in Chapter 3. Chapter 4 predicts the global dissipation of tidal energy at M2 due to 
seamounts. Conclusions follow in Chapter 5. 
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Chapter 2 



Data Acquisition 



The travel-time differences of pulse-like signals are measured at 3 minute intervals in 
1988 and 2 minute intervals in 1989 between the source at Oahu and the four receivers 
(Figure 1-1). Travel- times are intermittently measured for a total of between four and 
eight months (Figure 2-1). The transmission distances are about 4000 km and the transit 
times are about 40 minutes. 

Spiesberger et al. [1989a, b] describes the maximum likelihood estimator of phase 
differences used to estimate differences in travel-time. The basic process is outlined in 
Appendix D. With this estimator, the change in acoustic phase at 133 Hz, 6<f>, is related 
to the travel-time difference, St, by, 



27 v6(f> 
C U>A t ’ 



where u> c = 27r • 133 (rad-s -1 ) is the center radian frequency of the transmitted signal, u> is 
the radian frequency of the travel-time oscillation, and At is the geophysical time interval 
over which travel- time differences are measured (184 s in 1988, 122 s in 1989) [Spiesberger 
et al., 1989a]. 
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1988 1989 1990 



Figure 2-1: Shaded regions indicate the times when tomographic data were recorded at 
each of the receivers. 
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Chapter 3 



Data Analysis 



3.1 Barotropic Model for Acoustic travel-time 

The change in acoustic travel-time due to barotropic tides is given by, 

AT M (<,r B ) = ^-[/ n&t-'^-ds + h f )<k|, (3.i) 

c 0 Wr R c o ^Tr C 0 J 

where the geophysical time at the receiver is t, the geodesic between the source and 

receiver is Tr, and the mean sound speed is c Q . The barotropic current and surface 

displacement at location s on the geodesic at time t are denoted by u(s, t — and 

respectively. The incremental on the geodesic is ds , and l(s) is the distance 

from the receiver. The constant is ki = 3.9 X 10~ 3 [Spiesberger et al., 1989b]. The effects 

of finite propagation time are taken into account by the correction term (~^)- The finite 

propagation delay amounts to a phase correction of order ±10° at semidiurnal periods. 

The first term in (3.1) contains the effect of propagating sound with and against 
the tidal current. The second term, which amounts to a 1 or 2 ms correction, contains 
the combined effects of vertical iso-therm displacement and pressure fluctuations which 
accompany the rise and fall of the sea surface. For barotropic tides the travel- time change 
is about the same for different acoustic multipaths [Spiesberger et al., 1989a]. 

The elevation, £, is obtained from the Naval Surface Weapons Center (NSWC) tide 
model derived by Schwiderski [1978]. Unfortunately, except for A/2, NSWC models of the 
tidal current field are unavailable. We therefore use a linearized version of Schwiderski’s 
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[1978] Laplace tidal equations (LTE) to predict tidal currents. The equations are, 



W = McosOv + [^(ou? - /3C)] - k d u, 



(3.2) 



§7 = - 2 ft cos $u - 3C)] - k d v, 

where A and 6 are longitude and colatitude respectively, u and v are east and north 
current speeds respectively, rj is Newton’s equilibrium tide (which is proportional to the 
astronomical tide generating potential), g is the earth’s gravitational attraction, R e is 
the radius of a spherical earth, SI is the Earth’s angular rate of rotation, and kj is a 
damping coefficient. The constants a and (3 derived by Schwiderski [1978] correct for the 
solid-earth tide, ocean tidal loading of the solid earth, and the associated gravity potential 
perturbations of all tides and loadings. The solution to (3.2) for u and v is outlined in 
Appendix A. 

Equation (3.1) is decomposed into contributions from each tidal frequency, u?*., of 
interest, and solutions are determined by collecting and integrating sin and cos 

terms separately. The resulting travel-time change model is, 

K 

AT ht {t,T n ) = ^ e k (Tn)cosu k t + f fc (r R )sinu^, (3.3) 

k = l 



where e*.(rR.) and ^(Tr) are the solutions to the separated integrals, and K is the total 
number of tidal frequencies calculated. 



3.2 Barotropic Model Comparisons with Observed Tidal 
Signals 

3.2.1 Tidal Signal Estimation 

Power-spectra of the differences in acoustic-phase exhibit prominent variability at tidal 
periods [Spiesberger et al., 1989a, b] (Figure 3-1). For each receiver, tidal signals are 
estimated by fitting the 8 most prominent tidal components to the data (Appendix A, 
Table A.l). Specifically, we determine a/, and b^; k = 1 , 2 , . . . , 8 , in, 

8 

H(t n ,T R ) - ^^ki^cosuktn -\-b k (Tn)su\u k t n ; n = 1,2, . . . , A, (3.4) 
k = 1 
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Figure 3-1: Power spectral density estimate for the 133 Hz phase-differences at receiver 
R2. This periodogram uses a Bartlett data window with 42 days of data obtained in 1988. 
The tidal peaks at the K 1, Af2, and S 2 frequencies have rms values which correspond to 
travel-time change amplitudes of 8.6, 23.5, and 12.3 milliseconds respectively. 

by linear least squares [Lawson and Hanson , 1974], where the phase-differences are mea- 
sured at times t n , and the total number of measurements equals N. Depending on the 
receiver, N equals 55,000 to 135,000, so the least squares problem is overdetermined. 
Each least-squares fit to phase-differences is compared to barotropic model travel-times 
to determine how much of the signal can be attributed to barotropic tidal effects. 

3.2.2 Comparison with Model 

Assuming all of the (3.4) signal is due to barotropic influences, each tidal component 
is proportional to the time derivative of the corresponding (3.3) model. Comparisons 
between observations and the barotropic model are illustrated by the phasors of Figure 
3-2. Each phasor magnitude corresponds to a data-to-model signal amplitude ratio, which 
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Frequency 


Rl 

Amp(ms) 


R2 

Amp(ms) 


R3 

Amp(ms) 


R4 

Amp(ms) 


M2 


26 


22 


15 


7 


S2 


11 


7 


4 


3 


Kl 


10 


7 


1 


3 



Table 3.1: Predicted acoustic signal amplitudes for barotropic tides. The model ampli- 
tudes, (e£ + f^) 1 / 2 , are determined from the constants in (3.3) and have rms errors of 
about 2 ms. 



from (2.1) is, 



2fr( a t + b %) 1/2 



u c iv k A t{e\ + fl) 1 / 2 ’ 

and each phasor angle corresponds to a difference between data and model signal phases, 
which is given by, 



tan' 



■'(©-“-(sM- 



3.2.3 Analysis 

First, we discuss the comparison at two semidiurnal periods (M2 and S2). Predicted tidal 
amplitudes vary from 3 to 26 ms at our receivers (Table 3.1). At receivers Rl, R2, and 
R3, observed amplitudes are 0.8 to 1.6 times predicted amplitudes (Figure 3-2). At R4 
the observed amplitudes are 2.4 and 3.4 times predictions. Predicted and observed phases 
are within ±30° at Rl, R2, and R3 but differ by up to 47° at R4. Data at Rl, R2, and R3 
agree with predictions from the barotropic model to first order, but at R4, the barotropic 
model is inadequate. We return to this problem in the next section. 

We next discuss the comparisons for one semidiurnal period (K 1). The models for 
barotropic currents at diurnal periods are inaccurate within a ±10° band centered at 30° 
N (Appendix A). The transmissions to all receivers intersect this band, so comparisons 
between predictions and observations are problematic. Nonetheless, predictions at Rl 
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Figure 3-2: Comparison of phase and amplitude of measured and predicted tidal signals 
at four receivers (R1 to R4) and three frequencies (M2, S 2, K 1). Predicted tidal signals 
are obtained from the barotropic model given by (3.1). Each vector angle corresponds to 
a measurement phase that leads the predicted phase by the indicated amount. The unit 
circle in the upper left hand corner illustrates the vector angle convention. Vector magni- 
tude corresponds to the ratio of measured amplitude to predicted amplitude. Predicted 
magnitudes are listed in Table 3.1. An exact match between measurement and prediction 
is represented by a vector of angle zero and length one (tip just touching the unit circle at 
zero degrees). Solid arrows correspond to measurements of the combined 1988 and 1989 
record, and the dashed arrows correspond to separate 1988 and 1989 data measurements. 
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and R.2 resemble observations at tidal period K1 (Figure 3-2). It may be that R1 and 
R2 tidal signals are similar to predictions because their transmission paths travel smaller 
distances within 30° ± 10° N than paths to R3 and R4 (Figure 1). 

3.3 Baroclinic Tide Prediction 

3.3.1 Phase-Stable Internal Tides 

Apart from generation regions, baroclinic tides are not phase-locked to the ATGF, and 
their phases are incoherent at 100 km distance (approximately an internal tide wavelength) 
[Hendershott, 1981]. In generation regions, baroclinic tides may be phase-locked to the 
ATGF [Munk et al., 1981]. Open ocean evidence of this phase-lock is seen in a phase- 
stable M2 baroclinic tide over a large guyot east of Hawaii [Noble et al., 1988]. The 
transmission path to R4 passes within a few hundred kilometers of several large seamounts, 
one of which is indicated in Figure 1-1. We hypothesize that the semidiurnal signals at 
R4 do not agree with predictions from the barotropic model because of the presence of 
phase-locked baroclinic tides near the seamounts. 

The transmission paths to Rl, R2, and R3 may be near other large seamounts, but 
because there are no large discrepancies between the data and barotropic model predic- 
tions, it is unlikely that internal tides play a zero order role in determining tidal signals 
at these receivers. 

Appendix B summarizes a model for the generation of baroclinic tides from the inter- 
action between barotropic currents and seamounts. For baroclinic tides, acoustic travel- 
times are modulated by sound-speed fluctuations associated with vertical straining of 
iso-tachs [Spiesberger et al., 1989b]. The effect for the jth acoustic multipath, Fr(j), is, 

AT^TrU)) = i / Cint (r,i - (3.5) 

c o Jr K (j) c 0 oz 

where the position along the multipath is r, the vertical displacement of the internal tide 
is Cint? and the vertical derivative of the potential sound-speed is ♦ 



18 





M2 

Amplitude(ms) 


Phase(°) 


S 2 

Amplitude(ms) 


Phase(°) 


Barotropic Model 


7 


113 


3 


324 


Baroclinic Model 


20 


117 


5 


265 


Combined Model 


27 


116 


7 


286 


Measured Data 


24 ±1 


114 ± 1 


7 ±0.1 


279 ± 1 



Table 3.2: Comparison of predicted and observed travel-time changes at M 2 and S2 for 
receiver R4. The uncertainties on the measured data are based on the differences between 
1988 and 1989 least-square fits. Phases are based on a sin ut time dependence with t — 0 
at 0000 Greenwich mean time on 1 January 1989. 

3.3.2 Influence on Travel-Time at Receiver R4 

Being the southern-most receiver, R4 has a more east-west oriented ray path than the other 
receivers. This orientation results in a weaker influence from the generally north-south 
oriented M2 and S 2 barotropic tidal currents, so the effects of phase-locked baroclinic 
tides are more apparent. The predicted M2 baroclinic tide, generated from one of several 
seamounts near the transmission path to R4 (Figure 1-1), has an amplitude and wave- 
length of 6.2 m and 98 km respectively (Appendix B) and produces a 21 ms amplitude 
travel-time signal (Appendix C). 

The sum of the predicted M2 baroclinic tidal signals, generated from all prominent 
seamounts along the transmission path, is added to the barotropic prediction. The 
barotropic plus baroclinic prediction for R4 is almost identical to observations at M2 
and S 2 (Figure 3-3 and Table 3.2). The prediction is improved at K2 but not at N 2 
(Table 3.3). The frequency at N2 is close to M2 (Table A.l), and the observed amplitude 
at N 2 is small. We suspect that the barotropic plus baroclinic prediction does not fit 
observations at N 2, because strong baroclinic tides at M2 may be smeared into N 2 by 
dynamic processes [Hendershott, 1981]. 
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Barotropic Combined 

Model Model 



M2 



( \ 





Figure 3-3: The phasors on the left are identical to the combined 1988 and 1989 R4 record 
phasors shown in Figure 3-2. The phasors on the right present the same analysis using 
the combined barotropic and baroclinic models in Table 3.2. 





N2 

Amplitude(ms) 


Phase(°) 


K2 

Amplitude(ms) 


Phase(°) 


Barotropic Model 


1 


37 


1 


153 


Baroclinic Model 


6 


315 


3 


117 


Combined Model 


6 


325 


4 


126 


Measured Data 


2±1 


331 ± 38 


5±2 


153 ± 18 



Table 3.3: Same as Table 3.2, except comparisons are at N 2 and K2 frequencies. 
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Chapter 4 



Tidal Dissipation 



The effect of the ocean tidal bulge on the Earth-Moon system is a well known phenomena 
that is discussed extensively by Lambeck [1980]. The friction-induced phase lag between 
£ and 77 produces a torque in the Earth- Moon system that is responsible for the length- 
ening of the day and recession of the Moon at rates of about 1.8 ms-century -1 and 3.7 
m-century -1 respectively [Christodoulidis, et al. 1988]. The applied torque transfers an 
amount of orbital momentum to the Moon equal to the rotational momentum lost by the 
Earth [ Lambeck , 1977]. The net result of this momentum transfer is an overall reduction 
in system kinetic and potential energy. System energy is dissipated for the most part 
in oceanic tides, with some fraction (which has yet to be quantified) occurring through 
conversion to baroclinic tides [Lambeck, 1980]. In this chapter we estimate the barotropic- 
to-baroclinic conversion rate of the M 2 tide at seamounts. 



4.1 Tidal Dissipation at M2 from Satellite Observations 

The rate of energy dissipation is, 

-§■ = " 9 - 58 ^) 10 “ Oer g • s_ ‘> ( 4 -!) 

where is the total tidal acceleration of the Earth and 4j± \ s the acceleration of the 
Moon’s orbit in rad • s -2 [ Lambeck , 1977]. Recent estimates for M2 accelerations have 
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been obtained from the orbital perturbations of 17 satellites, and have values given by, 

— -4.45 x 10 -22 rad • s -2 , 

dt (4.2) 

% = -9.73 X 10“ 24 rad ■ s -2 , 

[Christodoulidis, et al., 1988]. Substituting (4.2) into (4.1) yields 2.5 X 10 19 erg-s -1 at M2 
(about 75 % of all tidal dissipation). 

4.2 Tidal Dissipation by Seamounts 

We consider two models for the distributions of seamounts. In the first model the number 
of seamounts with a basal radius r > R (km) is given by, 

N (R) = 1.44 x 10 6 exp[— 0.63/Z]. (4.3) 

Equation (4.3) is an extrapolation of Jordan et al. ’s [1983] results, for the eastern Pacific, 
to the world ocean. In the second model, the number of seamounts is estimated from 
altimeter measurements by Seasat [ Craig and Sandwell , 1988]. Seasat detects many large 
seamounts which are absent in the model by Jordan et al. [1983], but it does not detect 
small seamounts, as the altimeter profile sensitivity to seamounts falls off rapidly below a 
basal radius of 20 km [Craig et al ., 1988]. These two seamount distributions are outlined 
in Table 4.1. 

The rate of energy conversion from barotropic to baroclinic tides at each seamount is 
approximately, 

dF 

— = 2(A x)c gi Ei, (4.4) 

where Ax is the horizontal dimension of each seamount perpendicular to the direction of 
barotropic current, c g i is the group speed of the internal tide (Equation (B.22)), and Ei is 
the average energy density of the internal tide (Equation (B.21)). The dimension, Ax, is 
obtained from the seamount’s basal radius according to Table 4.2. If we take the barotropic 
current speed to equal 2 cm-s' 1 and the inertial frequency to equal 0.73 X 10“ 4 s _1 (30° 
N), then equations (4.4) and (B.23) yield a total energy conversion rate of 2 X 10 16 and 
1 X 10 18 erg-s -1 for the exponential and Seasat distributions respectively. These estimates 
are proportional to Ax and the squares of both seamount height and the barotropic 
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Basal Radius 
r(km) 


Number of Seamounts 


Exponential 

Model 


Seasat 

Model 


1-2 


360, 000 


0 


2-4 


290,000 


0 


4-6 


83,000 


0 


6-10 


30,000 


200 


10-15 


2500 


1200 


15-20 


110 


6000 


20-26 


5 


8800 


26-32 


0 


6800 


32-38 


0 


4800 


38-44 


0 


2800 


44-50 


0 


1600 



Table 4.1: Seamount distribution for the world’s oceans as a function of basal radius by 
Jordan et al. [1983] (exponential model) and Craig and Sandwell [1988] (Seasat model). 
The Seasat model is obtained by converting Craig and Sandwell’s [1988] histogram of 
peak- to- trough separation, p, to basal radius using r=(0.83)p in accordance with mea- 
sured seamount slopes. The numbers so obtained are multiplied by four to account for 
undetected seamounts [Craig et al., 1988]. Some of the discrepancy between models may 
be due to large seamounts in the western Pacific that were not represented in the Jordan 
et al. [1983] survey. 



Rectangular Seamount Dimensions 





Exponential 


Seasat 


Dimension 


Model 


Model 


Ax 


(0.62)r 


(0.62)r 


Az 


(0.21)r 


2500 m 



Table 4.2: Seamount dimensions for the exponential and Seasat models (Table 4.1). The 
length and width of each seamount is Ax , and the height is Az. Basal radius is denoted 
by r. The scale factor for Ax , (0.62), is derived from a mean flatness of 0.31 reported for 
85 seamounts in the Pacific [Smith, 1988]. The Az scale factor for the exponential model 
is derived from the mean height-to-basal-radius ratio of 0.21 [Smith, 1988]. The height is 
set to 2500 m for the Seasat model. 
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current speed (Equation (B.23)). The exponential model yields insignificant dissipation 
because it contains relatively few wide seamounts (Table 4.1). The Seasat model yields a 
dissipation that is about 4 % of the total at M2. Munk [1966] and Bell [1975] obtained 
rates of 5.0 X 10 18 and 2.5 X 10 18 erg-s -1 respectively for the conversion to baroclinic tides 
by bottom roughness. These estimates are the same order as we obtain using the Seasat 
distribution. 
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Chapter 5 



Conclusions 



We reach the following conclusions: 

1. Because barotropic signal predictions at Rl, R2, and R3 compare favorably with 
observations and because barotropic wave-lengths have basin-scales, the principle 
discrepancy between barotropic signal predictions and observations at R4 is probably 
not due to errors in Schwiderski’s [1978] model. Furthermore, because barotropic 
plus baroclinic signal predictions compare favorably with observations at R4, we 
believe the Schwiderski [1978] model is adequate for barotropic predictions. 

2. Seamounts may be responsible for a small (~ 5 %), but non-negligible, fraction of 
the dissipation of barotropic tides. The dissipation is sensitive to the heights, radii, 
and numbers of seamounts (all poorly known). 

3. The generation of baroclinic tides at seamounts can be accounted for with a two- 
layer model of the ocean (Appendix B). 

4. Barotropic tidal currents can be predicted from Schwiderski’s [1978] tidal elevation 
predictions with our model (Appendix A), except within ±10° of 30° North and 
South (diurnal frequencies) and near 90° N (semidiurnal frequencies). 

5. Comparison of measured tidal signals (from basin-scale acoustic transmissions) with 
predicted tidal signals (from barotropic models) can be used to detect baroclinic 
tides which are phase-locked to the ATGF. 
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6. It is difficult to make a-priori predictions of tidal signals in basin-scale transmissions, 
because generation regions of baroclinic tides are poorly mapped. 
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Appendix A 



Obtaining Tidal Currents from 
Elevations 



We solve (3.2) for u and v by defining, 



u = b\ cos ut + 62 sin u)t, 
v = a\ cos ut + a 2 sine 



(A.l) 



and, 

Trib [if 01 ? - w] = * cosut + » sina,( ’ (A 2) 

-ft [jj(a»7 - /?0] = *1 cos ^ + z 2 sina;i, 

where a; is the tidal frequency of interest and a \ , a 2 , &i» b 2 -> 3 / 1 , 2/2 1 Zi, and z 2 are constants. 
The elevation, £, and Newton’s equilibrium tide, 77 , are periodic with frequencies shown 
in Table A.l. The y’s and z’s are computed from (A.2) by expanding known functions, ( 
and 77 , in terms of sines and cosines. Substituting (A.l) and (A.2) in (3.2) leads to four 
linear equations, in the unknowns a\, u 2 , 61 , and & 2 , which can be written in matrix form 
as, 



0 —20, cos 9 —uj 

—uj kd 0 

— 20 cos 9 0 kd 



k t \(a l \ 



20 cos 9 

UJ 



&2 

&1 



\ kd 



uj 20 cos 9 



0 



\ b 2 J 



f V2 \ 
2 2 
yi 

\ z i ) 



(A. 3) 
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Semidiurnal Modes 


Diurnal Modes 


Component 


Frequency 

(cpd) 


Component 


Frequency 

(cpd) 


M2 


1.93227 


K1 


1.00273 


S2 


2.00000 


01 


0.92954 


N2 


1.89600 


PI 


0.99726 


K2 


2.00547 


Q1 


0.89325 



Table A.l: Major Tidal Frequencies 



The solution to this matrix equation determines u and v. 

We use the values, 

a = .69, 0 = .90, (A.4) 

derived by Schwiderski [1978], and the value of the linear damping coefficient, k di is 
selected by minimizing the phase difference between predictions of AT(/, M2, F/a) using 
Schwiderski [1983] and damped LTE currents. The resulting best value is, 



k d = 9.0 x 10 -6 s -1 . 



(A.5) 



The value of k d also controls the Q of the LTE model. We estimate Q by comparing 
the equations of motion to a forced linear-damped oscillator, 



d 2 x dx 2 

-Qp + k— + u 0 x = A cos wi, 



(A.6) 



whose Q is defined as, 




where the resonant frequency is, 



ur = {<4 - 




(A.7) 



(A.8) 



[ Marion , 1970]. By letting east be the x-direction, and substituting (A. 2) into (3.2), we 
obtain the following x-direction momentum equation, 
d 2 x Ox 

-Q ~2 + kd- 7 £ ~ fv = 2/1 cosut + j/ 2 sinu?i, (A.9) 
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where / = 2ficos 0 is the inertial frequency. Since, 



|f | ss 1.4 X 10 -6 , 

\fv\ « 1.4 X 10" 6 , 

{y\ + y 2 ) 1 ^ 2 ~ 4.4 x 10 -7 , 

\k d u\ « 9.0 x 10~ 8 , 

the dominant terms in (A.9) can be written 

du 

m~ fv = 0 - 

Substituting the solution to (A. 11), 



lull 




and the relation, 



u = icox, 



into (A9), we obtain the approximation, 



d 2 x dx 

Ht* + ^ i ~dt + lj2x ~ yi cosu:t + 3ft sino;t. 



(A. 10) 



(A. 11) 



(A.12) 



(A.13) 



(A.14) 



Comparing (A.6,7,8) with (A.14) we see. 






f ) l/2 



k d 



u _ 1.41 x lO'^- 1 
T d ~ 9.0 X 10 -6 s -1 



= 16, 



(A. 15) 



where 1.41 X 10 -4 s -1 is the frequency of the M 2 tide. This agrees with estimates ranging 
from 5 to 25 provided by Hendershott [1981]. 

A comparison of the damped LTE model with the Schwiderski [1983] model for M2 
in the North Pacific ocean demonstrates good agreement (Figure A-l). Both models also 
generally agree with the data presented by Luyten and Stommel [1990] extracted from 
deep moored current meter records. By extension, we believe that the damped LTE 
model should be an accurate first order predictor of the S 2, N 2, and K2 semidiurnal 
currents. 

The damped LTE model does not perform as well as the Schwiderski [1983] M2 current 
model in the vicinity of amphidromic points. The finite difference approximations we use 
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Figure A-l: The M 2 barotropic tidal current in the northeast Pacific Ocean at 0600 
Greenwich mean time on 1 January 1989 as determined from (a) the Schwiclerski [1983] 
model (solid arrows) (b) linear-damped Laplace tidal equations (Equation (3.2); dashed 
arrows) (c) data extracted from deep moored current meter records (two arrows desig- 
nated by circles around their tails). The current meter analysis is taken from Luyten and 
Stommel [1990] 
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for the spatial derivatives in (3.2) are not accurate near these areas of rapidly changing 
phase. These inaccuracies are limited to a 100 km (the grid size of Schwiderski’s model) 
radius, and the net effect on predicted travel-time change is less than 5 % in amplitude 
and 3 degrees in phase for a 4000 km geodesic passing through an amphidrome. 

The damped LTE model has an additional difficulty for diurnal currents within a ±10° 
band centered at the inertial latitude, 

6j = sin _1 (cjfc/2fi). (A. 16) 



For diurnal tides, Uk ~ 1 cpd and Of = 30° N, and for semidiurnal tides, Uk ~ 2 cpd and 
Of = 90° N. The difficulty can be seen by setting k d = 0 and solving (3.2) for u and v. 
The resulting equations, 



5 [-2ficose^(a7j-/3C)+^ ! |^(Q'^-/30] 

.R(4Q 2 cos 2 6 — u> 2 ) 

(A. 1 7 ) 



.R(4f2 2 cos 2 $— a/ 2 ) ’ 

have denominators equal to zero at the inertial latitude. These zeros are canceled by 
numerator zeros when w, v , and ( are determined exactly from the same equations of mo- 
tion, but a singularity occurs when the numerators are approximated by finite differences. 
The finite kj we use in (3.2) removes the singularity in (A. 17), but the solution is still 
inaccurate near 30° N (Figure A-2). 
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Figure A-2: The K 1 barotropic tidal current in the northeast Pacific Ocean at 0600 
Greenwich mean time on 1 January 1989 as determined from the linear-damped Laplace 
tidal equations (Equation (A. 2)). Note the irregular performance of the model in the 
vicinity of 30.1° N. 
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Appendix B 



Guyot- Generated Internal Tides 



In order to determine the coupling between barotropic tidal currents and internal tides 
generated by guyots, we use the two-dimensional two-layer model illustrated by Figure B- 
1. This model is applicable to sea plateaus whose major horizontal extent is perpendicular 
to the direction of the barotropic current. We follow the procedure used in deriving 
the coupling between surface tides and coastal generated internal tides [ Rattray , I960]. 
Variables are a function only of the distance from the center of the plateau (y-direction) 
in this two-dimensional model. 

A rigid-lid approximation is used. The momentum equations are, 



d t v' + fu f = 0, 



(B.l) 



d,v" + fu" = -g r d y <;", 



(B.2) 



and, 



d t u' - fv' = 0, 
d t u" - fv" = 0, 



(B.3) 



where u is the current parallel to the plateau, (') terms correspond to upper layer and (") 
to lower layer terms, / is the inertial frequency, and, 



9r = 9&P/P, 



(B.4) 



is reduced gravity, where, 



A p = P" - /»'. 



(B.5) 
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Figure B-l: Definition Figure for a two-layer ocean model used to predict the coupling 
between a barotropic tidal current and the internal tide generated by a seamount whose 
extent perpendicular to the direction of the current is considered infinite. 



Subtracting (B.2) from (B.l), we find, 

d t (v' - v") + f(u’ - u") = g r d y (", (B.6) 



and the continuity equations for each layer are, 



dy{h'v') = dtC , 
dy{h"vT') = -d t ("- 



(B.7) 



Let v s = (a/d) exp i(ut + <f ) ), which is essentially constant for any distance y small com- 
pared to the wavelength of the surface tide. We have contributions to the currents in the 
two layers from both the barotropic current and the currents associated with the internal 
tide. Thus, 



v' = v s + v[ u' = u s + u[ 
v" = v s + v” u" — u s + u", 



and from the rigid-lid approximation, 



(B.8) 



C = o, 

C" = Ct 



(B.9) 
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We choose the following solutions which satisfy (B.3) and (B.7), 



y<-( 

Ct -Be 1 *" 

V s l/d 2 

v\ -(Lj/k 2 h')Be ik *y 
v'/ (u/k 2 h")Be ik ™ 
u s —iffud 2 
u\ (if/k 2 h')Be ik *y 
u'l -(if/k 2 h%)Be ik *y 



-i < y < i 
A sin k\y 

i M 

— ( i(jj/kih')A cos k\y 
(iw/kih")Acos k\y 
-if/udx 

—(f/k\h')A cos k\y 
{f /k 2 h") A cos k\y 



V > t 
Be~ ik ™ 
l/d 2 

-(cj/k 2 h')Be- ik *y 
(. u/k 2 h%)Be- ik *y 
—iffu)d 2 
(if/k 2 h')Be- ik *y 
-( if/k 2 h")Be - ik ^ , 



with the common multiplier, 

a exp*^' 1 '^ . 



(B.10) 



(B-l 1) 



To satisfy (B.6), 



, 2 _ uA-B ' • _ i o 

S “ g r h'h'] a 3' 3 ~ L >*' 



(13. 12) 



In order to patch the solutions at the edges of the plateau, (" must be continuous at y = £ 
and y = —£. Thus, 



Be~ ik * 1 = A sin k t l. 



(B.13) 



The solutions (B.10) already satisfy the conservation of total mass flux h'v' + h"v n = 
dv s + h'v'i + h"v" across y = i and y — —t, but conservation of p' mass flux requires that 
h'v' = h'(v s + v t -) be continuous also. To satisfy this requirement we let, 



A = Aoe~’ (o “' /2) , 

B = £ 0 e- i(o,_,r/ 2-W). 



(B.14) 



Substituting (B.14) into (B.13) leads to, 



Aq sin k\t = Bq. 



(B.15) 



Using these relations and conservation of p' mass flux we find, 
h'cj . , h' iu 

— — ^4o(sin» + i cos a) sin k\t = — — j4o(sina -f icosa)coski£. (B.16) 

d 2 k 2 d\ k\ 

By equating the imaginary part of (B.16) we obtain, 

tan a = — tan k\£, (B.17) 

k 2 
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M <e 






\y\>( 


Amplitude 


Phase(°) 


Amplitude 


Phase(°) 


PC 

'sint 


aAo sin k\y 


0 — a + 90 


a Bo(\y\/y) 


<£-a + 90-fc 2 (|y|-^) 


V s 


ajd\ 


0 


a/d 2 


4> 


< 


aw Ao cos kiy/ki h! 


0 — a 


awBo/k2h l 


<)> - a + 270 - k 2 (\y\ - i) 


v" 


owAq cos kiy/kih" 


0 — a + 180 


awBo/k 2 h2 


</>-<* + 90 -fc 2 (|y|-f) 


U s 


af/wdi 


<£ + 270 


af/wd 2 


<t>+ 270 


u'i 


afAo cos kiy/ki h! 


0-0 + 270 


afBo/k 2 h ' 


</.-«+ 180 -k 2 (\y\-l) 


u" 


afAo cos kiy/kih'J 


0 — a + 90 


afBo/hh'J 


<t>-oc-k 2 {\y\-l) 



Table B.l: Generalized two-layer internal wave solution for a barotropic north-south tidal 
current propagating over an undersea plateau of infinite east-west extent. 



and by equating the real part, we can write, 



1 1 kih'cosa 

Ao = (- — ) j— 7, 

d 2 «i w cos kit 



(B.18) 



which leads to, 

Bq = (~ - -j-)^-sino. (B.19) 

d 2 d\ W 

The result of this solution is a standing internal wave over the extent of the plateau which 
patches to traveling internal waves at either end. The generalized solution is collected in 
Table B.l, and the particular solution for aguyot located at 30.3°N, 222.5° E is summarized 
in Table B.2. 

The cycle-average baroclinic energy density for this two-layer ocean is, 



Ti = i £ l p[iCmt)2gr + ft ' W)2 + ft ' W)2 + h " {u : ' )2 + ( B - 2 °) 

where T is the period of the internal tide. Substituting the Table B.l amplitudes and 
phases in (B.20) results in, 



Ei = ^pg r (v s d 2 B 0 ) 2 



(B.21) 



- P 

where the resulting energy density is in joules per square meter. This baroclinic energy 
propagates at a group speed that from (B.12) is given by, 

iP-p 



k 2 w 



(B.22) 
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Constituient 


M2 


S2 


Amplitude 


Greenwich 

Phase(°) 


Amplitude 


Greenwich 

Phase(°) 


A" 

( iint 


6.21 m 


166 + k 2 y 


2.08 m 


202 + k 2 y 


Vs 


2.3 cm-s -1 


68 


0.7 cm-s -1 


104 


v[ 


1.9 cm-s -1 


346 + k 2 y 


0.6 cm-s -1 


22 + k 2 y 


v” 


0.3 cm-s -1 


166 + k 2 y 


0.1 cm-s -1 


202 + k 2 y 


U s 


1.2 cm-s -1 


158 


0.4 cm-s -1 


194 


U 'x 


1.0 cm-s -1 


76 + k 2 y 


0.3 cm-s -1 


112 + k 2 y 


< 


0.2 cm-s -1 


256 + k 2 y 


0.1 cm-s -1 


292 + k 2 y 



Table B.2: Two-layer M2 and S2 internal tide solutions for regions north of a 22 by 111 
km undersea plateau. Seamount is located at 30.3°N, 222. 5°E (Figure 1-1), t — 11 km, 
di = 2200 m, d 2 = 4700 m, h' = 700 m, g T = .0059 m-s -2 , k 2 ( M2) = 6.41xl0 -5 m -1 
(wavelength is 98 km), and k 2 (S2) = 6.71xl0 -5 m -1 (wavelength is 94 km). 



The group speed, energy density, and horizontal extent of each seamount are combined 
in (4.4) to obtain the barotropic-to-baroclinic tidal energy conversion rate. Expansion of 
(4.4) in terms of fundamental parameters results in, 



where, 




(B.23) 



(B.24) 



The internal tide in Table B.2 has an energy density and group speed of about 1.5 X 10 5 
erg • cm -2 and 160 cm-s -1 respectively. The energy density near this seamount is about 
17 % of that contained in the parent barotropic tide, and the baroclinic conversion rate 
is about 5.3 X 10 14 erg-s -1 (53 Mw). 

The reduced gravity (B.4) used in the dispersion relation (B.12) is subject to the 
choices for the representative densities p' and p" in the upper and lower layers. Values of 
g r of order 10 -2 m-s -2 are typical for two-layer ocean models, but using g T = 10 -2 m-s -2 
is not optimal for minimizing the mean square error between the combined model and the 
measured travel time signals at R4. The optimal value of reduced gravity chosen for this 
model, as illustrated by Figure B2, is about 5.9 X 10 -3 m-s -2 . 
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Reduced Gravity (m/s 2 ) 

Figure B-2: Magnitude of Receiver 4 barotropic plus baroclinic signal prediction errors. 
The errors are plotted as a function of reduced gravity. 
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Appendix C 



Baroclinic Tide Travel-Time 
Model 



Equation (3.5) predicts the travel-time modulation of a single multipath caused by one 
seamount. The multipaths oscillate about the sound axis, and a typical wavelength (upper 
plus lower loop) is 50 km. The potential sound speed gradient in the region of the Table 
B.2 plateau is approximately, 



dcp ~ 
dz ~ 



0.13 s 1 above axis 
0.0 s -1 below axis, 



(C.l) 



where positive 2 is upward. We assume the average ray path is above axis for 50 % of its 
travel adjacent to the 111 km length of the plateau, so the average change in travel- time 
at the M2 period is, 



AT bc (t) = 



(6.2e Mt 4) )+ *<int m )(0.13s -1 )(lll km) 
(2)(1479 m s -1 ) 2 



: 21 e 






(C.2) 



where 0^ int is the phase of the M2 internal tide, and u is the M2 frequency. 

An evaluation of N s seamount-generated tides along the geodesic between source and 
receiver will produce N s baroclinic signals at each tidal frequency. The vector sum of 
baroclinic models plus the barotropic model results in the combined model, 



AT tot (t,u> k ) - AT ht (t,oj k ) -f Y2 A7bc(n,*,Wfc), 



(C.3) 
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where ATb c (rc, t,Uk) is the baroclinic signal prediction for the nth seamount at frequency 



^k- 
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Appendix D 



Signal Processing 



D.l The Transmitted Signal 



Pulse compression techniques described by Spiesberger et al. [1980] and Metzger [1983] are 
utilized to increase signal-to-noise ratio without degrading the phase difference resolution. 
The transmitted signal has a center frequency of 133 Hz and a bandwidth of about 17 Hz. 
The center frequency is phase modulated every 8 cycles using a 511 digit maximal shift 
register sequence having a period of 511(8/133) ~ 30.7 s [Spiesberger et al., 1989a]. The 
ideal transmitted waveform is, 



[Spiesberger et al., 1980], where f c is the center frequency, r ^ represents the shift register 
sequence of ones and zeros, and, 




(D.l) 




(D.2) 



where u(Z) is the unit step function. 
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D.2 Post-Detection Processing 

D.2.1 Demodulation and Averaging 

The received signal, z(t ), can be approximated by, 

N mp 

z(t) = Y, M‘ - Tn), (D.3) 

n=l 

which is the weighted sum of N mp multipaths each having gain b n and delay r n [Spiesberger 
et al., 1980]. Demodulation of the received signal is achieved by sampling at 4 • 133 Hz = 
532 Hz. Samples of an ideal single-path arrival are, 

z[n] = cos( Y) cos (|(r*[ n] - 1/2)] + sin( ) sin[^-(r/t[ n ] - 1/2)], (D.4) 

where, 

oo 

*[n] = X)u[n-32/], (D.5) 

/=0 

and u[n] is the unit step sequence. Because the sampling frequency is four times the 
carrier frequency, the even numbered samples are the cosine of the demodulated phase 
angle (weighted by ±1) and the odd numbered samples are the sine of the demodulated 
phase angle (weighted by ±1). The demodulated phase angles are low-pass filtered by 
block averaging to achieve a sequence of complex demodulates defined as, 

c[j] = x[j\ -f iy[j]\ j = 0, 1, . . . , J - 1, (D.6) 

where the square root of minus one equals z, 

*[?]= Y (— l) m+ M2(m + 4j)], (D.7) 

m =0 

and, 

y[j] = ^(-l) m+ ^[2(m + 4j)+ 1], (D.8) 

771=0 

[Spiesberger et al., 1980]. The resulting sample rate is four complex demodulates per 
digit, so there are J = 4 • 511 = 2044 complex demodulates per sequence. The complex 
demodulates are averaged over 6 consecutive sequences in 1988 resulting in one averaged 
sequence record every 184 seconds. The demodulates are averaged over 4 sequences in 
1989. The signal-to-noise ratio increases by 101og 10 6 = 8 dB following the six- period 
average in 1988 and by 101og 10 4 = 6 dB following the four-period average in 1989. 
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D.2.2 Replica Correlation 

One period of a 66.5 Hz sampling of (D.l) phase angles (expressed as complex demodu- 
lates) is contained in the following sequence of length J, 

= c °s[^(rjt|j] - i)] + tsin[|(r t[i j - i)]; j = 0,1,..., J - 1, (D.9) 

where, 

%'] = X>b'-4'h (D.10) 

/=0 

and u [j] is the unit step sequence. The received complex demodulates (D.6) are circularly 
convolved with a sequence-reversed (matched filter) replica of the transmitted demodu- 
lates (D.9). The pulse-compressed output is given by, 

J - 1 

d[m] = X] c[j']sJ[j + rn]; m = 0, 1, . . . , J - 1, (DTI) 

3=0 

where the periodic sequence, s^[j], is given by, 

oo 

»4>\j\ = H + (D.12) 

r=— oo 

[Oppenheim and Schafer, 1989], and s J is the complex conjugate of s The replica 
correlation (D.l 1 ) compresses 30.7 s of energy along each multipath into a single pulse 
of one digit length (~ 0.06 s) [Spiesberger et al., 1989a], whose amplitude and location 
depend on the gain, 6 n , and delay, r n , in (D.3). The signal-to-noise ratio increases by 
10l°gi O 511 = 27 dB following pulse compression, and a sliding average of 4 complex 
samples (separated by 0.015 s) increases the signal-to-noise ratio by another 6 dB at the 
cost of some resolution [Spiesberger et al., 1989a]. 

D.2.3 Phase Change Between Sequence Records 

The change in phase of the mth pulse-compressed demodulate from record r — 1 to r is 
given by the phase of, 

d[m, r — 1 \d*[m, r]. 

The average change in phase from record r — 1 to r is given by the phase of 

M 

A e tS< f>[ r ] = ^ d[m,r — l]cP[m,r], (D.13) 

m= 0 
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where there are M demodulates containing a signal [Spiesberger et ai, 1989a]. The phase 
in (D.13) is referred to as the cross-correlated phase, because it is the phase of the complex 
cross correlation of adjacent pulse-compressed records at zero-lag. The cross-correlated 
phase can be converted to the average travel-time difference between records using (2.1). 
The precision of the resulting estimation of travel-time difference is about 135 ps [Spies- 
berger et ai, 1989a]. 
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